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Abstract. It is well known that the surface brightness of elliptical galaxies and of bulges of spiral galaxies is best 



fitted by the Sersic 



profile. It is thus interesting to explore the lensing properties of the Sersic model because 



of its wide range of applicability. To this aim, we evaluate the lensing potential, the deflection angle, the time 
delay between the images and the amplification for a circularly symmetric lens whose surface density is described 
by the Sersic profile. We estimate the same quantities also adding a shear term to the defiection potential in 
order to study (at the lowest order) the impact of deviations from radial symmetry or the contribution from other 
nearby lenses. Moreover, we investigate the systematic errors due to the use of a de Vaucouleurs profile instead of 
the correct Sersic model also taking into account the presence of a dark halo. 
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1. Introduction 

Elliptical galaxies present a striking regularity in their 
global luminosity distributions in the sense that, within 
a wide range of sizes, their surface brightness pro- 
files can be described by simple fimctions. The de 
Vaucouleurs or r^/* law Ijde Vaucouleurs 1948|l was first 
proposed as a quite general function to fit the light pro- 
file of the elliptical galaxies and of the bulges of lentic- 
ular and spiral galaxies. However, it was soon realized 
that the Sersic r^^"^ profile, which generalizes the de 
Vaucouleurs law, is best suited to describe the surface 
brightness distribution of these systems IjUaon et al. 19931 
Kj^raham fc (Jolless 19971 [Prugniel fe Simien 1997| ). The 
regular properties of ellipticals have been the subject of 
different approaches, concerning both photometric and 
spectroscopic parameters, resulting in interesting scaling 
relations. Some well known examples are the fundamen- 
tal plane relating the effective radius rg, the luminos- 
ity intensity /g at rg and the central velocity dispersion 
fTo HDjorgovski fc Davis 1987| IDresseler et al. 19 871: the 
photometric plane linking the three parameters of the 
Sersic profile (^Khosroshahi et al. 20001 |Graham 2002|l and 
the entropic plane ( |Marquez et al. 2001| ). 

The regularity of elliptical galaxies is an attractive fea- 
ture for cosmological applications of gravitational lens- 
ing since it makes it possible to reduce the degener- 
acy among lens parameters and hence the systematic 
uncertainty in the estimate of cosmological quantities. 
The time delay among the images in multiply imaged 
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quasars systems may be used as a tool to determine the 
Hubble constant Hq IjRefsdal 1964|l once a model for the 
lens has been fitted to the observed images configura- 
tion, flux and time delay ratios. The use of scaling re- 
lations makes it possible to narrow the parameter space 
expressing some of them as function of the remaining 
ones. On the other hand, it has been estimated that at 
least 80% of gravitational lenses are elliptical galaxies 
HTurner et al. 1984L |F?ikugita fc Turner 1991J . Hence, it 
is important that statistical lensing studies, such as the 
number counts of lens systems, the distribution of image 
angular separations and of time delays, take into account 
the most correct model of elliptical galaxies in the anal- 
ysis of the available data. Up to now, this has not been 
done since these studies have been concentrated on the 
dark halos rather than on the luminous components of 
the galaxies. Therefore, this latter component is usually 
neglected or simply modeled with the spherical Hernquist 
profile ( [Hernquist 1990| ) which approximately reproduces 
the de Vaucouleurs law when projected to the lens plane. 
The only remarkable exception is the work by Maoz & 
Rix (1993) where the lens is modeled using the r^/'^law 
directly, together with an isothermal dark halo. It should 
be interesting to repeat their analysis using the more real- 
istic Sersic profile and the scaling relations we have quoted 
above. As a first step, one has to study the lensing prop- 
erties of this model. This is the aim of the current work. 

The paper is organized as follows. In Sect. 2, we deter- 
mine the lensing potential and the deflection angle for a 
circularly symmetric lens described by the Sersic model. 
We are thus able to write the lens equations whose solu- 
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tions are found in Sect. 3. Section 4 is devoted to the eval- 
uation of the amplification of the images, while the effect 
of a shear term, due to deviations from circular symmetry 
or to nearby lenses, is investigated in Sect. 5. A study of 
the systematic errors due to the use of a de Vaucouleurs 
profile mimicking the best - fit Sersic model to the lens sur- 
face brightness is presented in Sect. 6. We summarize and 
conclude in Sect. 7. 



respectively. Using the above change of variable, Eq.Q) 
becomes : 

dx^ ax 



in which : 



(6) 



2. Lensing potential and deflection angle 

The surface mass density of a galaxy whose surface bright- 
ness is fitted by the Sersic model may be written as 
IjSersic 1968|) : 



S(r) T /e exp <^ -b{n) 



l/n 



(1) 



in which T is the mass - to - light ratio, /g the luminos- 
ity density at the effective radius and h(n) a constant 
defined such that the luminosity within is half the 
total luminosity. It is possible to show that b{n) may 
be found by solving the following equation (see, e.g., 
IMazure fc Capelato 200"2| ) : 



r(2n,fe) = r(2n)/2 



(2) 



where r(a, z) is the incomplete F function and r(a) the 
actual r function. 

Adopting a rectangular coordinate system in the lens 
plane centred on the lens galaxy and with {xi,X2) ori- 
ented along its main axes, r is simply {x\ -\- ^2)^^^. We 
will assume that the surface density is circularly symmet- 
ric so that all the lensing quantities will depend only on 
r. While useful in the computations, this approximation 
is not a serious limitation to our analysis since the re- 
sults for the circular case may be immediately generalized 
to flattened models by means of numerical integration 
HSchramm 19901 IKeeton 200 Moreover, we will inves- 
tigate later the impact of deviations from circular sym- 
metry by adding a shear term to the lensing potential. 
However, we stress that the results for the circularly sym- 
metric models allow us to obtain a picture of the main 
properties of the Sersic model as a lens. 

It is convenient to introduce a new dimensionless vari- 
able defined as : 



l/r, 



(3) 



The lensing potential ipir) may be found by solving the 
following equation l|Schneider et al. 1992|l : 



(4) 



where K(r) — S/Ecrit is the convergence with Scrit = 
c^Ds/4:TtGDiDis the critical density for lensing and Dg, 
Di, Dis are the angular diameter distances between ob- 
server and source, observer and lens and lens and source, 



Eq. (|SJ) may be easily solved using the Mathematica pack- 
age. We obtains : 



^(x) = Kg— x2"2^2[{2n,2n},{l-f 2n,l + 2n},-6x] 

-n K^e^ b-^"rl T{1 + 2n) \nx + ci Inx + C2 , (7) 

where pFg[{ai, . . . , Oj,}, {61, . . . , 6^}, y) is the generalized 
hypergeometric function^ ( jGradshteyn fc Ryzhik 1980| ) 
and ci and C2 are two integration constants that have to 
be fixed. Since the lensing potential is defined modulo an 
additive constant with no physical meaning, we can put 
C2 = 0. As regards ci, we may impose the constraint that 
the potential is regular both at zero and at infinity, i.e. : 

^{x) ^ a; ^ 0, cxD . 

This condition allows us to choose the constant ci in such 
a way that the final expression for the lensing potential 
turns out to be : 



2„ 2-F2[{2n, 2n}, {1 + 2n,l + 2n}, -bx] 
2F2[{2n, 2n}, {1 + 2n, 1 + 2n}, -b] 



with the value of the potential for r = rg 
given by^ : 



-He 2^2[{2n,2n},{l + 2n, l + 2n},-6] . 



(8) 

X = 1, 

(9) 



Given the circular symmetry of the model, the deflection 
angle will be purely radial and has an amplitude a given 
by dip /dr. Using the new variable x and differentiating, 
we get : 



a{x) 



e X 



= 27ireKefo"^"e''a:~" r(2n) 



[r(l + 2n)-2nr(2n,6x)] 
r{2n,bx) 



1 - 



r(2n) 



having used the relation r(l + 2n) ~ 2nr{2n). Using 
Eq.Qj we may finally write the following expression for 
the deflection angle : 



a{x) = 2af. X 



1 



r(2n,&x) 



r(2n) 



(10) 



^ We use here the same notation for the generalized hyper- 
geometric function as in the Mathematica package. 

^ By using Eq.JUJ, it is possible to eliminate the denominator 
in Eq.Q. However, we prefer to write the deflection potential 
as in Eq.© so that the exponent n is the only model parameter 
entering the scaled deflection potential V'/V'e- 
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Fig. 1. The scaled deflection angle a/cke vs the dimen- 
sionless radius r/re for three values of the exponent n of 
the Sersic profile, i.e. n = 3 (long dashed), n = 4 (solid), 
n — b (short dashed). 




Fig. 2. The characteristic deflection angle ae vs the lens 
redshift zi for three values of the exponent n of the Sersic 
profile, i.e. n = 3 (long dashed), n — A (solid), n = 5 
(short dashed). 



in which : 

ae = a{r = re) = nre Keb~'^"e''T{2n) 



(11) 



Note that we are assuming that the distances in the lens 
plane are measured in arcsec so that a has the dimen- 
sion of arcsec and the lensing potential of arcsec^. Ea. fTIHl 
should be obtained also directly since, for a circularly sym- 
metric model, it is l|Schncider ct al. 1992.. Kceton 2001,1 : 



a{r) = — 



2 r S(r') , , , 
r ar 



^ Jo ^crit 

2nre KeX-^" / exp [~b{x' ~ l)]dx' 

Jo 

r{2n,bxy 



2ae X 



1 



r(2n) 



(12) 



which is the same as Ea. (|10|l . Note that the agreement be- 
tween the results for a from the two different procedures 
is an indirect check that our choice of the two integration 
constants ci and C2 in Eq.Q is indeed correct. It is inter- 
esting to observe that only the exponent n of the Sersic 
profile determines the behaviour with the dimensionless 
radius x of the scaled deflection angle a/ae^ while the ef- 
fective radius re, the luminosity density le at re and the 
mass - to - light ratio T enter only as scaling factors. As an 
example, Fig.^ shows a{x)/ae for three values of n. 

This plot shows two interesting qualitative features. 
First, we note that a/ae quickly decreases with r/re with 
the larger values of n going to zero faster. This is an ex- 
pected result since it simply reflects the fact that the sur- 
face mass density of the Sersic model is exponentially de- 
creasing with r^/". On the other hand, we also observe 
that the shape of a/ae is almost independent on n for 
values oi r/re > 1 as is apparent from Fig.^ We have 
checked that the plots for larger values of n are more and 
more similar as n increases, the only difference being the 
height of the peak. 



The scaling parameter Ue is a function of the Sersic 
model parameters n, re and /g and of the mass - to - light 
ratio T for a given cosmological model and lens {zi) and 
source (z^) redshifts. We observe that cte oc T as is ex- 
pected since the higher the lens mass, the higher the bend- 
ing of light in a given position in the lens plane. To qualita- 
tively study the dependence of Ue on the other model pa- 
rameters and the lens redshift, it is important to consider 
lens systems that are as realistic as possible. In Appendix 
A we give a detailed description of how we choose the 
Sersic parameters for the lens galaxy taking also into ac- 
count its redshift. To evaluate ae, we have to compute 
the critical density Scrit- To this aim, we need to fix 
the cosmological parameters and the lens and source red- 
shifts to estimate the angular diameter distances Ds, Di, 
Dis. We adopt the flat ACDM model with (17™, I^a, h) = 
(0.3,0.7,0.72), with h = i/o/lOO km s'^ Mpc'^ The 
source redshift Zg is fixed as kxzi in which k ~ 3.7 is 
the median value of Zs/zi among the observed lens sys- 
tems^. Fig.|21 shows the characteristic deflection angle ae 
as function of the lens redshift for three values of n and 
T = 10. The increase of ae with zi simply reflects the low- 
ering of the critical density because of the higher value of 
the distances ratio Dg/DiDis. Fig.|21also shows that the 
lower n is, the higher is ae- However, it is difficult to give a 
qualitative explanation for this behaviour. Ea. (|ll() seems 
to suggest ae oc n6^^"r(n). Actually, the dependence of 
ae on n is much more complicated since n also determines 
b through Eq.ljSJ, while re is linked to n and zi given the 
way we have chosen the Sersic parameters. These consid- 
erations prevent us from giving any (at least) qualitative 
interpretation of Fig.|3 



3. The lens equation 

The time delay of a light ray deflected by the galaxy lens- 
ing effect is given by : 



See the CASTLES web page IIKochanek et ai. 200211 . 
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Fig. 3. The Einstein radius vs the lens redshift zi for 
three values of the exponent n of the Sersic profile, i.e. 
n = 3 (long dashed), n = 4 (solid), n — h (short dashed). 



-r^ - rrs cos {9 - 9^) + -r^ - tp{r, ( 



(13) 



where {r,9) is the image position, {rs,9s) the unknown 
source position and ipir, 9) the lensing potential and we 
have defined : 



DiDA (1 



(14) 



According to the Fermat principle, the images lie at the 
minima of At, so that the lens equations may be simply 
obtained by minimizing At. Inserting Eg.© into Ea. lfT^ 
and differentiating, we get : 



r — Ts cos {9 — 9s) — 2aeX 
Ts sin (0-0,) = . 



1 



T{2n,bx) 



r(2n) 



= 0, (15) 
(16) 



Ea. ljlGI) has two immediate solutions. The first is — 0, 
i.e. the source and the lens are perfectly aligned. In this 
case, the image is just the Einstein ring with radius 
obtained solving Ea. lfT^ for — 0. We plot in Fig.|31 
the Einstein radius as a function of the lens redshift 
for three values of n and T = 10. Note that the te 
depends significantly on n for lenses with redshift ap- 
proximately in the range (0.2,0.6). This is just the red- 
shift range of the most of the observed lens systems (see, 
e.g., IKochanek et al. 2002|l thus suggesting that the use 
of the correct value of n is quite important, i.e. using the 
de Vaucouleurs profile instead of the Sersic model may lead 
to serious systematic errors. We turn back to this problem 
in Sect. 6. 

The second solution is obtained by imposing 
sm{9-9s) = 0, i.e. 9 = 9^ + rmr (with m = 0,1). In 
this case, we get two images symmetrically disposed with 
respect to the lens centre. The radial coordinate ri of the 





Fig. 4. The radial coordinates ri and r2 of the first (top 
panel) and second (bottom panel) image vs th source po- 
sition Vs for three values of the exponent n of the Sersic 
profile, i.e. n = 3 (long dashed), n — 4 (solid), n — 5 
(short dashed), with the mass-to- hght ratio T — 10. 

first image (i.e., the one with 9 = 9s) is obtained by solv- 
ing Ea. ((TB|l with cos {9 — 9s) — 1, while the second has 
a distance r2 from the lens centre obtained by solving 
Eq.l^nil with cos {9 - 9s) = -1. In Fig.gl we plot ri and 
r2 as function of the source radial position for the different 
values of n used. The Sersic parameters and the cosmolog- 
ical background are fixed as usual, while we have adopted 
izi,Zs) = (0.31,1.7 2), as f or the real lens PG 1115-^080 
dWeymann et al., 1 980 Im pey et al., 1996| ). 

Note that we get two images which is an even number, 
while one could invoke the theorem on the odd number 
of images IjSchneider et al. 19921 IStraumann 1998|l to sus- 
pect that (at least) a third image should appear. However, 
the theorem is not violated since it only holds if the mass 
distribution and its first derivative have no singularities. 
This is not the case for the Sersic model. While S(r) in 
Eq.(^ is regular, its first derivative is singular in the origin 
and hence the above theorem does not apply. 

An important observable in lens systems is the time 
delay between the two images of a distant source. This 
may be simply evaluated as : 



Ati2 = At(ri,0i)-At(r2,02) 



(17) 



V. F. Cardone: The lensing properties of the Sersic model 



5 



Dti2 HdaysL 



0.05 0.075 


0.125 


0.15 


0.175 


0.2 


^^^^ « 











Fig. 5. Time delay between the two images of a distant 
source as a function of the source position for three 
values of the exponent n of the Sersic profile, i.e. n = 3 
(long dashed), n — i (solid), n — 5 (short dashed). 
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Using Ea. lfT^ . we may quickly compute the time delay 
between the two images formed by a circularly symmet- 
ric Sersic lens. This is plotted in Fig.|Slas function of the 
source radial coordinate for a system with {zi,Zs) = 
(0.31,1.72) and T = 10. Note that the time delay scales 
with the mass - to - light ratio as since it is proportional 
to the square of the image positions and the latter (as the 
Einstein radius) scales linearly with T. Fig. shows that 
the first image always arrives first being less delayed by 
the bending of light due to the lens. On the other hand, 
the time delay is quite short when compared to the val- 
ues usually measured, the lowest being ^ 10 d. However, 
drawing any conclusion by this comparison is not possi- 
ble since one should also take care of the observational 
problems related to the measure of short time delays. 

4. The amplification of the images 

The bending of light due to gravitational lensing changes 
the solid angle subtended by a source on the sky and leads 
to an amplification of the images with respect to the unde- 
fiected case. The image amplification is simply the inverse 
of the Jacobian of the lens mapping. It is : 



A{xi,X2) 



1 



detJ 



= [(l-V^ll)(l-V^22)-^?5 



(18) 



in which 7/;^ = d^Tp/dxidxj. For a circularly symmetric 
lens, this reduces to (ISchneider et al. 1992|l : 



A{r) 



1 '^^11 '^'^ 
r J \ dr 



1 - 



r^x" 



1 - ^2;!-" — 
nrp dx 



1 -1 



(19) 



Inserting Ea. (|10l) into Ea. (|19|l . with some simple algebra, 
we finally get : 



A{x) - [{l-p){l+p + q)Y 



(20) 



7 



].iai 



Fig. 6. The amplification of the first (top panel) and sec- 
ond (bottom panel) image as a function of the source dis- 
tance from the lens centre for n = 3 (long dashed), n = 4 
(sohd), n = 5 (short dashed). 



having defined : 



p{x) 



q{x) = 



2ae 



2a. 



r(2n,6a;) 
r(2n) 



nre r(2n) 



exp {—bx) 



(21) 



(22) 



The amplification of the two images is shown in Fig.|Bl 
where the Sersic model parameters, the cosmological back- 
ground and the lens and source redshifts have been fixed 
as for Fig.^ 

This plot shows that, while the first image is always 
amplified and conserves the same orientation as the source 
(i.e., it has positive parity), the second image always has 
negative parity and could also be deamplified (i.e., the 
magnification is lower than 1). Therefore, we conclude that 
for certain values of n and rs a circular Sersic model is 
unable to produce two observable images and hence its 
lensing effect is actually null. 

It is interesting to study the critical curves of this 
model. These are defined as the loci of the points in 
the lens plane such that the amplification becomes (for- 
mally) infinite, while their projection on the source plane 
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gives the caustics. Ea. (|20|) shows that A diverges either if 
l—p = or if 1+p+q = 0. The first condition individuates 
the tangential critical curves which, the model being cir- 
cularly symmetric, will be a circle with radius given by the 
solution of the equation 1 — p = 0. Inserting Eq.|2U into 
this condition, it is easy to show that the only tangential 
critical curve is just the Einstein ring. The corresponding 
caustic is thus the point = 0. On the other hand, the 
condition l+p + q = determines the radius of the radial 
critical curve. However, using Eas. (|21|l and (|22|l . this rela- 
tion reduces to an equation which has no solutions. Thus, 
we conclude that the circular Sersic model has no radial 
critical curve. 

5. The effect of the shear 

We have considered the case of a circularly symmetric 
model, but we are aware that real elliptical galaxies are of- 
ten flattened systems'*. If the deviations from circular sym- 
metry are not too large (i.e., if the eccentricity e = l — b/a 
is not too different from 1), the effect of flattening may be 
taken into account (to the lowest order) adding a shear 
term to the lensing potential which is now : 



^ = V'ic„s(r-) - ^-fr^ cos {29 - 26^) 



(23) 



in which '!/'icns('') is given by Eq.(|SJl, while 7 is the shear 
intensity and 9^ its orientation angle. A shear term could 
also mimic the contribution to the lensing potential of 
nearby lenses or of the cluster (if any) to which the main 
lens galaxy belongs. Adding the shear term to the lens- 
ing potential obviously changes the lens equations. Some 
simple algebra allows us to finally obtain : 



sm 







7 sin {20 - 29^) ' 

[r{9)-Q{e)y^+n{9)^o 

having defined : 

r{9) = [1 + 7 cos (261 - 261^)] sin^ {0 - 9,) , 
Q{9) = 7 cos {0-0,) sin [0-0,) sin (261 - 29^) 



n{9) 



n 2-2 

laeTe^ sm 



(261 - 261^ 



T{2n,hx) 
T{2n) 



(24) 
(25) 

(26) 
(27) 

(28) 



To find the images position, we can thus solve Ea. H25|l 
with respect to 9 for given values of the source coordinates 
{rs,9s) and of the shear parameters (7,^7) and then use 
Ea. H24l) to find the corresponding value of the radial co- 
ordinate, excluding those solutions of Eq. 125() which give 
rise to negative r. Without loss of generality, hereinafter, 
we will assume 0^ — 0. Moreover, we fix 7 = 0.11 as 

* Note, however, that bulges of spiral galaxies are typically 
more rounder than elliptical galaxies and hence the circular 
case we have studied is not completely unrealistic. 



Fig. 7. Zero -level contour plot of the left-hand side of 
Ea. (|25ll for a source on axis, i.e. 0, = 0, as a function of 
the radial source coordinate (in arcsec) and the image 
angular coordinate (in rad) excluding the region giving 
rise to negative r. Results are shown for models with n — 3 
(long dashed), n = 4 (solid), n = 5 (short dashed). 



determined by Holder & Schechter (2003) through a set 
of numerical simulations^. The Sersic model parameters, 
the background cosmological model and the redshifts of 
lens and source are fixed as in Sects. 3 and 4. Because of 
the breaking of radial symmetry due to the shear term, 
the solutions of the lens equations now depend not only 
on the radial, but also on the angular coordinate of the 
source. In Fig.[71we show the contour plot at the zero level 
of the left - hand side of Eq. (|25|l as a function of and 
assuming the source is on - axis, i.e. 0, = 0, and excluding 
the regions which give rise to negative r. For a given rs, 
the solution of Ea. H25l) may be graphically found as the 
intersection points of the vertical line = const with the 
plotted curves. Two images are formed reducing to only 
one (the Einstein ring) if — (not shown in the plot). 
When the source is off axis, the number of images is two 
or four according to the source position, as may be argued 
looking at Fig. |H1 which refers to the case 0s = tt/3. 

The presence of shear also changes the amplification as 
is easy to understand observing that, because of the bro- 
ken circular symmetry, A will now be function of both the 
polar coordinates {r,0). Actually, Ea. lfT^ does not hold 
anymore and we have to use Ea. (|18|l to estimate A{xi, X2)- 
As a consequence, the critical curves are changed too. 
Fig-El shows that now we get two critical curves. The ex- 



^ The value 7 = 0.11 refers to real quadruply imaged systems 
IHolder fc Schechter 20031 . while here we do not know in ad- 
vance what is the number of images. The systematic error we 
are making using this value does not affect our main results. 
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Fig. 8. Same as Fig.Qbut for Os = tt/3. Note that, for 
graphics - related problems, the sharp edges in the plot are 
not real the actual ones being more smoothed. Note that 
the contours (in particular the lower ones) for different 
values of n are almost superimposed. 
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Fig. 9. Critical curves in the lens plane {xi,X2) for the 
Sersic model modified by the addition of a shear term. The 
plots refer to n = 2 (top left), 3 (top right), 4 (bottom left), 
5 (bottom right). Distances in the lens plane are measured 
in arcsec. 



ternal curve is the tangential critical curve which is now 
deformed into an ellipse instead of a circle as in the case 
without shear. Moreover, a second critical curve^ appears 
inside the first corresponding to the radial curve which 
was absent in the circularly symmetric case. 

6. Systematic errors 

Even if the Sersic profile is well known to be 
a better fit to the surface brightness of ellipti- 
cal galaxies IjCaon et al. 19931 IGraham &: CoUess 19971 
[Prugniel fc Simien 1997| ), it is customary in lensing stud- 
ies to use the de Vaucouleurs r^/** law to model the visible 
component of the lens. It is thus important to investi- 
gate the systematic errors introduced by this procedure. 
To this end, we proceed as follows. First, we fix the Sersic 
parameters as described in the Appendix arbitrarily fixing 
the lens redshift zi = 0.31. As a second step, we generate 
a simulated set of surface brightness measurements using 
realistic error bars. These data are then fitted with the 
de Vaucouleurs profile to obtain the values of and le 
that best reproduce the observed surface brightness. This 
procedure ensures that the deVacouleurs parameters are 
reasonably chosen and the projected mass distribution of 
the lens galaxy is well represented both by the true Sersic 
profile and the fitted r^/'* law. 

® Note that the inner critical curve in the upper left plot in 
Fig.|^is rotated with respect to the other plots. It is important 
to stress that this is only an artifact of the numerical solution. 



As a first application, we consider the estimate of the 
mass - to - light ratio. It is quite easy to show that : 

j,obs 

^ = —( f T im (29) 

rE[n,re,Ie, T = 10) 

with r'^'^ the observed value of the Einstein ring^. Note 
that we have fixed T = 10 in the denominator in Ea. (|29|l 
just to make the numerical computation easier. Let us 
suppose now that we have seen an Einstein ring and we 
have been able to measure the lens surface brightness. 
Then, we can fit this profile using both the Sersic and 
de Vaucouleurs law so that we know the three parameters 
(n, re,/e) and may evaluate T using Ea. H29|l . Using the 
de Vaucouleurs law with n = 4 instead of the correct Sersic 
profile leads to a different estimate of T. Fig.^| shows 
log (Tdv/'^'sersic) as a function of the exponent n of the 
Sersic model used to generate the simulated data, being 
TTdv (Tsersic) the value obtained modeling the lens with 
the de Vaucouleurs (Sersic) law. We only include in the 
plot those models that are well fitted by both profiles. 

The results are quite interesting showing that the 
mass - to - light ratio is underestimated if n < 4 or over- 
estimated if n > 4 by a factor which may be as high as 10 
(excluding values of n < 2 or n > 5 that are less realis- 
tic). However, Fig.llOlshould not be considered as a defini- 
tive proof against the de Vaucouleurs model instead of the 

^ Actually, it is not so easy to observe an Einstein ring, but 
we may qualitatively estimate rs from the image separation in 
multiply imaged quasars. 



8 



V. F. Cardone: The lensing properties of the Sersic model 



logHUdV UsersicL 




n 



Fig. 10. The systematic error in the estimate of the mass- 
to - Hght ratio T as a function of the exponent n of the 
correct Sersic model describing the lens. 



Sersic one. Actually, this result have been obtained us- 
ing simulated observations of lens galaxies whose param- 
eters have been chosen with a procedure that could also 
be revised. It should be interesting to repeat this analysis 
for real lens systems also comparing the estimated mass- 
to- light ratios with those predicted by models of stellar 
evolution. 

In the previous discussion, we have assumed that the 
lens redshift is known which is not always the case. Fig.O 
shows that the way ve scales with zi depends on n too. 
Henceforth, we should correctly write : 



dV 



- Sersic 



(n = 4,re,/e,T = 10) 



„obs 



rf=-«'=(n,re,/e,T = 10) ' 



with rj^ (rlf^"^) the Einstein radius evaluated for the 
de Vaucouleurs (Sersic) model. Combining the previous 
two relations, we get : 



Tdv ^ r^'^'-^'^(n,re,/e,T = 10) 
Tscrsic r|V(n = 4,re,/e,T= 10) 



(30) 



This depends only on n if we know zi (as assumed above), 
otherwise it is a function of both z; and n. This suggests 
that the results in Fig. EH should be read with some cau- 
tion. Let us consider, for instance, the case rt = 3. From 
Fig.EIIwe get Tdv/TTsersic < 1. On the other hand. Fig. 13 
also shows that rlf ^"^ < rj^ for values of the lens redshift 
lower than about 0.2. As a result, we conclude that, if z/ 
were not known, we should be unable to say whether the 
systematic error is due to having chosen an incorrect lens 
model or a wrong redshift. Actually, the situation is still 
more involved since the way we have fixed the Sersic pa- 
rameters also implictely depends on zi . We thus stress the 
need to repeat this analysis using real galaxies to thor- 
oughly investigate the systematics induced by using the 
de Vaucouleurs profile instead of the Sersic one. This is 



outside the scope of this paper, but will be the subject of 
future work. 

Real lensing galaxies are usually modeled adding a 
dark halo to the visible component. If the Einstein radius 
has been measured somehow (e.g, by directly observing 
an Einstein ring or from the images separation), then it 
is possible to get an estimate of the projected mass of the 
dark halo inside the Einstein ring as : 



M^te) = 7rr| Y.crit - Mvis(f£;) 



(31) 



with Mvis (?■_£;) the projected mass of the visible component 
estimated from the Sersic (or de Vaucouleurs) model as : 



JQ 



2n-l^-b{x-l) 



dx 



(32) 



with xe = {te I'l'eY^"' ■ We may suppose that {n,re,Ie) 
have been measured fitting the Sersic law to the lens sur- 
face brightness profile and that T has been estimated from 
a stellar evolution synthesis code so that Eq.lPU may be 
used to directly obtain the projected mass of the dark 
halo. 

To investigate the systematic errors induced by the use 
of a de Vaucouleurs profile instead of the correct Sersic 
one, we first have to choose a dark halo model. We adopt 
the softened power law sphere (see 'Kee ton 200 11 and ref- 
erences therein) whose projected surface density is : 



S(r) 



2 (s2+r2)i-''/2 



= So 1 



2\ /3/2-1 



(33) 



which represents a flat core with scale radius s and cen- 
tral surface density Eqj and then a power law decline 
with exponent (3 defined such that the mass grows as 
asymptotically. The core radius can be zero if /3 > 0. 
As illustrative cases, we consider the two models with 
/3 = 1 (i.e., a softened isothermal model) and /3 = —2 
(i.e., a Plummer model). To fix the halo parameters, 
we first deproject both the halo and the Sersic pro- 
file (Binne y fc Tremaine 19871 [Mazure fc Capelato 2002D , 
then estimate the total mass of the visible component and 
fix the halo virial mass, Mvir, so that 90% of the total 
galaxy mass is represented by the halo itself. The value 
of the core radius s is then obtained by scaling that of 
the Milky Way as s = smw (Mvir/Af¥W)i/3^ ^1^^^ 
straightforward to evaluate Sq as function of Afvir and s. 

We use Ea. (|31|l to estimate the projected mass of the 
dark halo using for the lens visible component the Sersic 
model or the corresponding fitted de Vaucouleurs model 
assuming T = 10 for both profiles. The Einstein ring is ob- 
tained by solving Eq. (|15() with = and adding the halo 
contribution to the lensing potential. In Figs .1111 and 1121 
we plot log (-/Vfdv/-Wscrsic) sls sl function of the exponent n 
of the correct Sersic model, in which Afdv (A/scrsic) is the 
value of the projected mass of the dark halo obtained after 
the visible component have been fit by the de Vaucouleurs 
(Sersic) profile. As is apparent from Fig. ^2 the halo pro- 
jected mass for the case (3 = 1 may be seriously overes- 
timated (underestimated) if n < 4 (n > 4). Considering 
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Fig. 11. The systematic error in the estimate of the pro- 
jected dark halo mass inside the Einstein radius te as a 
function of the exponent n of the correct Sersic profile 
using the softened power law model with /3 ~ 1 for the 
halo. 



only models with n in the range (3, 5) which is more re- 
alistic, the error could be up to a factor 10 so that it is 
worthwhile to investigate more carefully this topic using 
real lens systems. 

On the other hand, the error is lower if the halo is mod- 
eled as a Plummer sphere as is shown in Fig. El However, 
the Plummer model is less suited to describe galaxy haloes 
so that it is not worthwhile to spend further work on this 
topic. 

7. Conclusions 

Simple analytical calculations show that elliptical galax- 
ies form 80% of the gravitational lenses. Hence, it is very 
interesting to investigate their lensing properties taking 
into account models which are as realistic as possible. The 
Sersic r^/" profile is well known to fit the surface bright- 
ness of both elliptical galaxies and bulges of lenticular and 
spiral galaxies. In this paper we have studied its lensing 
properties assuming circular symmetry of the surface mass 
distribution. 

Under this hypothesis, we have derived analytical ex- 
pressions for the lensing potential, the deflection angle and 
the image amplification. Solving the lens equations, we 
have found that a circularly symmetric Sersic lens pro- 
duces two images of the source that are merged in the 
Einstein ring if = 0, i.e. the source and the lens centre 
are perfectly aligned. One of the two images may be deam- 
plified and hence, for some choice of n, a Sersic lens does 
not produce two observable images, an effect that must 
be taken into account in statistical studies such as lens 
number counts. No radial critical curve is formed by a cir- 
cularly symmetric Sersic lens, while the tangential critical 
curve is the Einstein ring. 

Adding a shear term to the lensing potential allows us 
a first order investigation of the impact of deviations from 
circular symmetry and to take into account the contribu- 
tion of nearby lenses. We have solved the lens equations 



Sers 



\ 
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Fig. 12. Same as Fig.HTl but with (3 
model. 



-2 for the halo 



in this case assuming, without loss of generality, that the 
shear is on -axis. Due to the breaking of radial symmetry, 
the number and positions of the images formed depend 
now both on and 9s and not only on as in the sym- 
metric case. We have found that an on axis source gives 
rise to two images, while two or four images may appear 
as the source moves away from the lens axis. The presence 
of the shear also changes the amplification of the images 
and the shape of the critical curves. The tangential criti- 
cal curve is deformed into an ellipse, while an inner radial 
critical curve appears. The size of both curves depends 
on the Sersic model parameters and decreases with the 
exponent n. 

Our results suggest that a Sersic lens can produce 
up to four images of a distant source if a shear term is 
added to the potential. This is in agreement with the 
recent result by Kochanek (2003) who was able to fit 
constant M/L models to two quadruple QSOs, namely 
RXJ0911-I-0551 (Bade et al. 1997VBurud et al. 1998) and 
PG1115-I-080 CWeymann et al., 1980 , Impcy et al., 1996| ), 
using the de Vaucouleurs profile, that is to say a Sersic lens 
with n — A. However, a direct comparison is not possible. 
To fit RXJ0911+0551, one has to use two lensing galaxies 
modeled with the de Vaucouleurs profile, taking also into 
account the presence of the cluster which the lenses belong 
to. A similar situation also holds for PGl 115-1-080 where 
a single de Vaucouleurs lens is unable to fit the images 
configuration and the contribution of a nearby cluster is 
strongly required to reproduce the observable quantities. 

The Sersic profile has been shown to be a better fit 
to the surface brightness distribution of elliptical lens- 
ing galaxies than the de Vaucouleurs law. Nonetheless, 
the de Vaucouleurs model is usually used in lensing stud- 
ies. We have thus investigated the systematic errors in- 
duced by this procedure. This analysis has shown that 
the mass - to - light ratio T may be seriously under - or 
overestimated (depending on n being lower or larger than 
4) when considering constant AI/L models. On the other 
hand, if a dark halo is added to the visible component 
of the lens galaxy, the use of a de Vaucouleurs model 
instead of the correct Sersic one leads to a wrong esti- 
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mate of the dark halo projected mass of the dark halo in- 
side the Einstein ring. These considerations could suggest 
that the use of de Vaucouleurs profile to model galaxies 
that are best fitted by the Sersic law must be avoided. 
Actually, these results should be treated with great cau- 
tion since they have been obtained using simulated galax- 
ies whose model parameters have been fixed with a pro- 
cedure that is by no means unique. In particular, as de- 
scribed in detail in the Appendix, we have used the pho- 
tometric plane HKhosroshahi et al. 2000IIGraham 2002|l to 
relate the Sersic parameters (n, re,Ie)- Actually, this rela- 
tion has been determined from galaxies belonging to the 
Virgo and Fornax clusters, while we have extended it to 
higher redshift galaxies the validity of which is still to 
be proven. Henceforth, a different prescription to fix the 
model parameters could be implemented and a detailed 
study of how this affects the results is needed. However, 
the procedure we have adopted leads to Sersic parameters 
describing galaxies which look quite reasonable so that 
we are confident that our main results are trustworthy. 
Actually, the best strategy is to resort to a careful study 
of real lens systems in order to check whether the effects 
we have reported are indeed true, or a consequence of how 
we have chosen the model parameters. 

The analytical expressions of the lensing potential 
and of the deflection angle here obtained for a circularly 
symmetric Sersic lens may be easily generalized to the 
flattened case following the procedure described, e.g., in 
Keeton (2001). The resulting potential, deflection angle, 
and amplification may then be used to try to fit the ob- 
served lens systems (both double and quadruple) in order 
to investigate if constant M/L models are viable solutions 
or whether a dark matter component has to be invoked to 
reproduce the observable lensing quantities (images posi- 
tions and flux ratios). Since the Sersic model parameters 
(n,re,/e) may be directly measured from detailed pho- 
tometry, the only unknown lens quantity is the mass - to - 
light ratio T and hence we can study eventual correla- 
tions among T and the photometric parameters. On the 
other hand, we can assume that lensing galaxies obey the 
same scaling relations (such as the photometric plane and 
the fundamental plane) as the low redshift systems. If we 
model the lens as the sum of a luminous Sersic component 
and a dark halo, we can reduce the degeneracy among 
model parameters excluding those solutions which lead to 
Sersic parameters not allowed by the quoted scaling rela- 
tions. Moreover, if the time delays among the images were 
measured, this approach should reduce the systematic er- 
ror in the determination of the Hubble constant through 
the time delay method IjRefsdal 1964|l . On the other hand, 
statistics of gravitational lensing, such as lenses number 
counts and distribution of images angular separation, are 
known to be powerful tools to discriminate among differ- 
ent cosmological models. The modeling of the lenses plays 
a key role in this kind of analysis and it is thus very im- 
portant to use realistic lens models. The Sersic lens we 
have studied may thus be a first step towards a reanaly- 
sis of the available data in order to better investigate if it 



is possible to reconcile them with constant M/L models 
and if the constraints on the cosmological parameters are 
changed by the use of more realistic models of elliptical 
galaxies. These questions will be addressed in a series of 
forthcoming papers. 
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ing suggested us to investigate the systematics discussed in 
Sect. 6 and for the comments that have helped to improve the 
presentation. 

Appendix A: Choosing the lens parameters 

Let us now explain the way we choose the Sersic param- 
eters in order to have a lens system at redshift zi which 
is as realistic as possible. To this aim, let us first consider 
the Milky Way. Its visible component may be modeled as 
the sum of a bulge with total mass Mbuigc and an expo- 
nential disk with scalelenght Rd and local surface density 
'S{Rq) = Sq, in which Rq ~ 8.5 kpc is the Sun distance 
to the galactic centre. We set : 

Mbuigc = 1.3 X lO^*' Mq , So = 54 Mq pc'^ . 

We now build a galaxy with total intrinsic luminosity 
given as : 



Mtot = -2.51og('^^MdiskT 



(A.l) 



with Mtot the galaxy total absolute magnitude and 
Afdisk = 27ri?^ E0 exp{Ro/Rd) the total mass of the 
Milky Way* disk. The apparent total magnitude of this 
galaxy at redshift zi is then : 



^ltotizl) = Mtot + 25 + 5 log 



DLizi) 



Mpc 



+ K{zi) 



(A.2) 



with K(z) the if - correction (estimated by interpo- 
lating the values reported for elliptical galaxies in 
iPoggianti 1997§ and Dl the luminosity distance that, for 
the flat ACDM model, is : 



DUz)=^ii+z) r 

tia Jo 



dz' 



^nK + a„^{l + z'f 



On the other hand, /itot is related to the three Sersic pa- 
rameters as follows : 



^tot = -2.5 log 27rnrg/, 



(A.3) 



* Having chosen the Milky Way as reference for the lens mass 
scaling is not completely justified since our Galaxy is a spi- 
ral one, while we are considering elliptical galaxies as lenses. 
However, this does not introduce any systematic error in the 
results since what we really need is a realistic value for the 
mass of a galaxy. To this aim, we could take as reference any 
galaxy provided that the mass of its different components are 
well measured. We have thus chosen the Milky Way as refer- 
ence only because its mass has been better determined. 



V. F. Cardone: The lensirig properties of the Sersic model 



11 



with Te expressed in arcsec. Let us now reduce the Sersic 
model parameters space using the photometric plane as 
found by Khosroshahi et al. (2000; see also lGraham 2 002^1 
by fitting the Sersic law to a set of near - infrared observa- 
tions of elliptical and SO galaxies : 



0.172 logre - 0.069 hq = logn - 1.18 , 



(A.4) 



in which re is in kpc and fio is the central surface bright- 
ness related to le as follows : 



L = dex 



Ho + 1.08575(n) 



2.5 



(A.5) 



with dex(a;) = 10^. We have now all we need to describe 
our procedure of choosing the Sersic parameters for a re- 
alistic lens system. 

1. Choose a value for the mass - to - light ratio T and use 
Ea. ljA.l|) to estimate the absolute magnitude of the 
galaxy. 

2. Compute the apparent total magnitude at redshift zi 
using Ea. ljA.2|l . 

3. For a fixed value of n, solve Eas. HA.3|l . HA.4I) and (|A.5p 
with respect to r^ to get the value of the effective radius 
in kpc and convert it to arcsec. 

4. Estimate Jp as : 



lQ-l^tat{zi)/2.5 



(A.6) 



This simple procedure ensures us that the lens models 
we are using to investigate quantitatively the main lens- 
ing properties of the Sersic profile are meaningful so that 
the results we have obtained may be considered reliable. 
However, it is worth to note that we have implicitly as- 
sumed that the photometric plane coefficients in Eq. HA.4p 
do not change with redshift. Actually, this assumption has 
still to be verified, but we are confident that our main re- 
sults are not seriously affected by revisions of the photo- 
metric plane coefficients. One might also resort to a differ- 
ent procedure, not using Ea. (|A.4|) . For instance, one could 
follow the first two steps of our procedure, then fix both n 
and le and finally evaluate Tg solving Ea. (jA.3|) . However, 
following this proscription, wc find that the higher the lens 
redshift, the larger re in complete disagreement with what 
we observe (and with common sense) . This is not the case 
for our prescription. A detailed and quantitative study of 
how the results are affected by the procedure adopted to 
fix the model parameters is nonetheless needed and will 
be performed in a forthcoming paper. 
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